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In contrast to semiconductor structures, the experimentally observed plasma resonances in 
graphene show an asymmetrical and rather broad linewidth. We show that this can be explained by 
the linear electron energy dispersion in this material and is related to the violation of the generalized 
Kohn theorem in graphene. 

I. INTRODUCTION 

Plasma oscillations in semiconductor two-dimensional (2D) electron systems have been studied since late 1960-ies^— . 
Plasmon related absorption resonances have been observed in different semicondutor structures (inversion layers 4 -, 
quantum wells*, wires^— , dotsiSr— , antidot o 12 ' 14 ^ 15 , ring o 16 ^ 17 ) and in the broad frequency range from microwaves 
(the frequency / ~ 10 GHz) up to far-infrared (/ ~ 1 THz). Both the frequency and the linewidth of plasma 
resonances were found in good agreement with theoretical predictions. 

In graphene 1 ^— , a new truly 2D material with very attractive physical properties 2 ^, the spectrum of plasma 
waves was first calculated in 2006-200 7 22 ' 23 , see also 2 ^— . However, the number of experimental studies of graphene 
plasmons is not very large, in spite of great interest to this material. Their observation by the electron energy los o 34 i 35 
and far-infrared transmissio n 36 ! 37 spectroscopy has been reported only in few papers, see also^. In all experimental 
studies the linewidth of plasmon resonances was found to be comparable with or even larger than their frequency, 
although graphene samples with the mobility of 200000 cm 2 /Vs are currently available. 

Recently, it was theoretically shown 3 ^ that the width of the classical cyclotron resonance (CR) should be substan- 
tially broader in graphene than in semiconductors due to the strongly non-parabolic (linear) energy dispersion of 
graphene electrons 

e p = ±v F \p\ (i) 

(here vf ~ 10 8 cm/s is the Fermi velocity). In fact, it was found that the spectrum ([T]) leads to a broad observed CR 
linewidth even in perfectly pure graphene, in the absence of any scattering or other energy-loss mechanisms. The goal 
of this paper is to study whether the linear electronic spectrum ([T]) may be also responsible for the broadening of plasma 
resonances in graphene. Using a non-perturbative kinetic approach we study the motion of massless quasiparticles 
([T]) in a finite-width stripe of graphene with a parabolic confinement potential and show that, indeed, the observed 
plasmon resonance line should be strongly asymmetric and have the width comparable with the plasmon frequency. 
This qualitatively agrees with observations of Refsi^ir— . 

II. PLASMA OSCILLATIONS IN A GRAPHENE STRIPE 

The plasma oscillations in a solid can be modeled as a collective motion of electrons in an external potential 
produced by a positively charged background. For example, in a three-dimensional metallic slab occupying the area 
|x| < W/2, see inset to Figure QJi, the positive background with the charge density +en v 8(W/2 — \x\) creates the 
electric potential 



_ ™n v W 2 / (2x/W) 2 , \2x/W\<l, 
^ X >~ 2 \4|a;|/W-l, \2x/W\ > 1. ' (Z) 

shown in Figure QJi. Being shifted from their equilibrium positions in the x-direction, electrons of the slab move in 
the potential ([2|) in accordance with the equation of motion 

mi = J^- = ~4Tm v e 2 x, \x\ < W/2. (3) 
ox 

Since the potential inside the slab is parabolic, Eq. ([3]) describes harmonic oscillations with the frequency ui 2 3 = 
Aim v e 2 jra which coincides with the frequency of three-dimensional plasmons. 
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FIG. 1: The confining potential energy — e<j)(x) due to the positive background in (a) a 3D slab and (b) a 2D stripe. Insets 
illustrate the geometry. Red dashed curve in (b) shows the parabolic approximation ©. 



In a 2D stripe |x| < W/2, z = (inset to Figure[TjD) the positive background with the charge density +en s 6(W/2 - 
|a;|)(5(z) creates the electric potential 
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which can be also approximated, with a good accuracy, by a parabolic form inside the stripe, 

-e(t){x,z = Q)KiKx 2 /2, \x\<W/2, 
Figure [TJa. Choosing the potential curvature K from the condition 



W/2 



dcj)(x, z = 0) 
dx 



Kx J dx = 0, 



so that the approximation (J5)) creates the same average force inside the stripe, one gets 



K = 161n2 



(4) 



(5) 



(0) 



The motion of all electrons in the stripe is then described by the harmonic oscillator equation of the type ([3]) with 
the squared oscillation frequency ojp 2 = K/m — 161n2n s e 2 /mW , which differs by only 10% from the numerically 
calculated^ plasma frequency in a Unite-width 2D stripe. 

Due to the linear energy dispersion (JlJ the dynamics of graphene quasi-particles substantially differs from that of 
electrons in semiconductors. In the absence of scattering their oscillations in the field of the positive background can 
be described by equations 



vfp/p, 



F = (-Kx)e x 



(7) 
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with the curvature K given by ([5]). The problem ([7]) should be solved with initial conditions r(0) = i"o and p(0) = 
Po = Po(cos cj>o, sin0o)i where po < pf, i.e. we assume that the momenta of all particles at t = lie inside the Fermi 
circle. 

The system of equations ([7]) has two integrals of motion, the momentum p y o and the total energy 



Ho = vf^JpI +pl + K—= v fPo + K-±. (8) 



The time dependence of the coordinate x(t) is then determined by the formula 



t = ±l , n °- Kx ' 2 ' 2 dx >, (9) 



l*o v F ^(H -Kx> 2 /2)*-v F p y0 

the relation p x (t) can be found from (|8]), and the dependence y(t) - from the equation for y in ([7]). Equation (|9]) 
describes oscillations of the x-coordinate of the particle between two turning points xi.2 = ±y/2(Ho — v F\Pya\)/ K . 
The oscillation period T = T(p y0 ,Ho) = T(p ,(j) ,Xo) is then given by the formula 

T = 2 f n Q -KxV2 



v F ^(H -Kxy2y-v F p 2 y0 



Taking the limit xo — ¥ we finally get the oscillation frequency fl = 2-k/T in the regime of very small oscillation 
amplitudes: 



si(p M = \ i FK ,A v (ii) 
y po A(4>o) 



where the function 



2 /o ry/l-\ain<t>\ i _ 2 

A(0 = ^ / . dx (12) 

Jo - .x 2 ) 2 - sin 2 <j) 

has the properties = A{-(j>) = A(tt + cj>), A(Q) = 2V2/n rj 0.900, A(ir/2) = 1. The inverse function 

which determines the frequency (| 1 lj) is shown in Figure [2j 

In contrast to semiconductors, the oscillation frequency ft of massless quasiparticles, Eq. (jTTj) . depends on the 
initial conditions. It weakly varies with <pQ (Figure [5]) but quite essentially depends on the absolute value of the initial 
momentum po. Particles with the largest (Fermi) momentum have the lowest oscillation frequency 
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FIG. 3: The shape of the plasmon absorption line in a perfectly clean graphene stripe. 



see Figure [3] Particles with smaller momenta oscillate with higher frequencies, therefore the observed 2D plasmon 
resonance line should have a broad asymmetric form with a flatter high-frequency and a steeper low-frequency side, 
Figure [3] The maximum of the absorption plasmon curve lies at the frequency 



which can be treated as the 2D plasmon frequency in a graphene stripe. The formula (|14|) gives correct dependencies 
on all physical parameters and a reasonable numerical value for the plasma frequency in a graphene stripe (cf. witb4^) . 

The shape of the plasmon absorption line (Figure is very similar to the one that was experimentally observed in 
Refsj 36 i 37 . Similar to the classical cyclotron resonance 3 -^, the 2D plasmon line broadening takes place in the collisionless 
limit, in the absence of any scattering in the graphene sample. Taking into account the finite scattering rate would 
lead to a somewhat smoother curve and to an additional broadening of the line. The collisionless broadening of the 
2D plasmon resonance line is a direct consequence of the linear energy dispersion (p} and is related with the violation 
of the generalized Kohn theorem^— in electronic systems with a non-parabolic energy dispersion. 



Our results for the frequency of the 2D plasmon resonance (|14|) are in very good agreement with those of Refs . 22 ' 23 . 
As for the linewidth of the resonance, they are substantially different: our calculations give a finite- width resonance 
line in the collisionless limit, while those of Refs i 22 ' 23 give a vanishing linewidth under the same conditions. We believe 
that the reason of this difference results from the non-analyticity of the graphene spectrum ([lj in the Dirac point 
p = 0. The solution of the equations of motion ([7]) is equivalent to a non-perturbative solution of the Boltzmann 
kinetic equation. In Refsi^ 2 ^ 3 . the linear response approach was applied, i.e. the solution of the kinetic equation was 
expanded, from the very beginning, in powers of the electric field and only the first-order term was taken into account. 
Mathematically, such a procedure suggests that all coefficients of the kinetic equation are analytical functions of p, 
which is the case in semiconductor structures but is not the case in graphene. In order to fully understand the graphene 
response to external electromagnetic fields a more detailed non-perturbative classical (Boltzmann) or quantum theory 
is thus required. 

To summarize, we have predicted that the 2D plasmon resonance line should be asymmetric and substantially 
broadened, even in perfectly clean graphene (without disorder) due to the linear energy dispersion of graphene 
electrons. This broadening is not the case in semiconductor systems with a parabolic energy dispersion and is related 
to the violation of the generalized Kohn theorem in graphene. 
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